Quantum-to-classical crossover of quasi-bound states in open quantum systems 
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In the semiclassical limit of open ballistic quantum systems, we demonstrate the emergence of instantaneous 
decay modes guided by classical escape faster than the Ehrenfest time. The decay time of the associated quasi- 
bound states is smaller than the classical time of flight. The remaining long-lived quasi-bound states obey 
random-matrix statistics, renormalized in compliance with the recently proposed fractal Weyl law for open 
systems [W. T. Lu, S. Sridhar, and M. Zworski, Phys. Rev. Lett. 91, 154101 (2003)]. We vaUdate our theory 
numerically for a model system, the open kicked rotator. 
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Open quantum systems are described by non- self- adjoint 
operators and hence have different spectral properties than 
closed ones. Instead of an orthogonal set of stationary bound 
states at real energies one encounters mutually non-orthogonal 
quasi-bound states i/jn with complex energies En, which de- 
cay exponentially in time with a uniform decay rate 7n/2 = 
— Im E'n . The complex energies define the poles of the scat- 
tering matrix, for which the random-matrix theory (RMT) 
provides a universal benchmark in the case of random wave 
dynamics |T] (for applications to the transport in electronic 
nanostructures see |2]). RMT has been derived for disordered 
systems Ij], while its status for ballistic dynamics is unsettled. 
In this paper we investigate ballistic systems that remain open 
in the classical limit — the decay mechanism then corresponds 
to ballistic escape of point particles into the asymptotic scat- 
tering region. (This can be realized in clean nanostructures 
and in driven atoms, provided the escape does not require tun- 
neling.) One goal is to uncover and quantify deviations from 
RMT. We focus on the case of classically chaotic systems (our 
methods are easily transferred to integrable or mixed dynam- 
ics). 

We show that in approaching the classical limit of these 
ballistic systems, the quasi-bound states separate into two dif- 
ferent classes. The classes are discriminated by the Ehren- 
fest time TEhr = ^ In A/', up to which the escape observes 
quantum-to-classical correspondence f^. (N is the number of 
ballistic escape channels and A is of the order of the Lyapunov 
exponent.) Phase space regions with classical escape faster 
than TEhr support instantaneous ballistic decay with rate 7^ 
larger than the inverse time of flight through the system. The 
quasi-bound states with instantaneous decay develop a large 
degree of degeneracy, which is accompanied by a drastic de- 
parture from orthogonality — the states are almost-linearly de- 
pendent. Non-orthogonality of states has been identified jsllal 
as the principal obstacle for the formulation of Weyl laws to 
rigorously count the states in open systems. We restore the ap- 
plicability of the Weyl law and estimate the relative fraction of 
instantaneous decay modes as 1 — exp(— TEhr /^d well), where 



Tdweii is the mean dwell time. The overwhelming majority 
of long-lived quasi-bound states associated to phase space re- 
gions with classical escape slower than TEhr are amenable for 
an effective RMT fr], in which the dependence of the effective 
matrix dimension on Planck's constant conforms with the re- 
cently discovered fractal Weyl laws in classically chaotic open 
systems |6]. 

The motivation for our work arises from the up- surging 
interest in the quantum-to-classical correspondence of open 
quantum systems, witnessed in particular over the past two 
years. That quantum-to-classical correspondence invalidates 
the assumption of wave-chaotic motion and implies non- 
universal corrections RMT was pointed out some while ago 
LSJ . However, only recent experiments and numerical investi- 
gations could access the very large systems needed in view of 
the just-logarithmic scaling of TEhr with the effective system 
size (here given by N). In stationary transport, the Ehren- 
fest time manifests itself in the systematic suppression of shot 
noise |9, 10, 11] and in the on-set of classical conductance 
fluctuations 1 12] . It also governs the impact of decoherence on 
dynamical systems coupled to heat-bath environments (open 
systems of a different kind than investigated here) |^3l [ij]. 
Spectral footprints of the Ehrenfest time have only been ad- 
dressed for closed systems, in terms of the form factor of 
the spectral correlation function |8, 15] and the proximity- 
induced excitationgap in a normal conductor next to a su- 
perconductor BldllTl] . 

For simplicity, we formulate our theory of quasi-bound 
states for quantum maps F with one degree of freedom, oper- 
ating on states (pin 3. Hilbert space of finite dimension M 1 1 8] . 
This dimension serves as the inverse of the effective Planck's 
constant, /leff = l/M. The evolution over one unit time step 
including loss is given by (p{t + 1) = QF(p{t), where Q is a 
projection operator of rank M — N which describes survival in 
the system. This projection operator introduces sub-unitarity 
into the time evolution, which is the equivalent to non-self- 
adjointness in the energy domain. The mean dwell time is 

Tdwell = M/N. 




FIG. 1: (Color online) (a) Regions of escape after one (red), two 
(green), three (blue), four (magenta) iterations in the open standard 
map (|9j with K = 7.5 and escape for x G (0.2, 0.4) (rdweii = 5). 
The other panels show Husimi representations of quasi-bound states 
for Hilbert space dimension M = 160. (b,c) Short-lived states with 
instantaneous ballistic decay (|/i| < 10~^), localized on the classical 
pre-image of the opening. (d,e) Long-lived states with random- wave 
characteristics, (f) Trapped long-lived state (only one state of this 
kind exists for M = 160; it has the smallest decay rate). 



For M ^ oo (i.e., /leff -^ 0) and fixed Tdweih the system 
attains its classical limit, in which it is described by an area- 
preserving map JT operating on a bounded two-dimensional 
phase space of normalized area 1, while survival is described 
by a projection operator Q onto the complement of the open- 
ings. For illustration see the phase space portrait of the open 
classical standard map in Fig.[lla) (details of this system are 
given below). 

Quasi-bound states are defined by the condition of quasi- 
stationarity 



QFQl/jn = /J^ni^n, l^n = exp{-iEn), 



(1) 



of the internal part of the wave function. Since QFQ is a 
sub-unitary operator, the eigenvalues have modulus |/in| < 1, 
(i.e., lie inside the unit circle in the complex plane), and the 
quasi-energies En have a negative imaginary part. The decay 
rate of a state ?/^n is given by 7^ = — 2 Im En. 

Our objective is to identify and count instantaneous decay 
modes by exploring the strict quantum-to-classical correspon- 
dence observed for escape times shorter than the Ehrenfest 
time. At the outset, note the set of A^ trivial short-lived states 
in the kernel of Q, all of which have eigenvalue fim = 
(m = 1, . . . , A"). They can be collected into the rows of an 
M X A'-dimensional matrix Pq, which fulfills Q = t — PqPq . 
The basic building blocks for the construction of a much 
larger number of nontrivial short-lived states are the con- 
nected phase-space regions At^i of escape after t iterations 
(t = 1, 2, 3, . . .), depicted in Fig.[lla). We denote the union of 
all regions with fixed thy At = Ui^/t^t,*, where It contains 
all applicable indices. Formally associating a region ^0 to the 
opening of the system, these regions partition phase space. 

Next, we introduce the regions Bt = Ui^jj^At^i, where Ij. 



restricts the index to areas larger than a Planck cell /leff- In 
this part of phase space the escape of initial wave packets cor- 
responds to the classical particle dynamics. Consequently, we 
call B = Ut<rEhr^t the region of quantum-to-classical cor- 
respondence. The maximal index TEhr arises because of the 
finite size of the Planck cell, and defines the Ehrenfest time of 
quantum ballistic escape. 

In order to construct states supported by the region B 
of quantum-to-classical correspondence, let us introduce the 
characteristic projection operators lil9ll 



Vt = 



D{Bt) 



dxdp\x{x,p)){x{x,p)l 



(2) 



where \x{x.,p)) are minimal-uncertainty wave packets local- 
ized at position x and momentum p, and D{-) denotes the 
characteristic function of a region. In the strict limit /leff -^ 0. 
the operator Vt represents the characteristic function of the re- 
gion At- For finite /leff, only the regions Bt are well-resolved 
by the wave packets. The operators Vt are defined with these 
smaller regions since this guarantees idempotency up to small 
corrections due to leakage. Hence, each Vt projects onto some 
subspace Ht of a dimension Mt = dim. Ht. This implies 
the representation Vt = PtP^, where the M-dimensional 
columns of P^ are mutually orthogonal, VtPf = ^tt'Pt- We 
also introduce the complementary projector V = 1 — Xlt ^t = 

T 

PP . 

Due to quantum-to-classical correspondence in the regions 
Bt, the semiclassical dynamics propagates states from sub- 
space Ht to Ht-i and finally to the opening, where they are 
destroyed by Q\ 

QFQHtCHt-i (t>l); (3a) 

QFQHt = {0} (t = 0,l). (3b) 

Equation (l3bl immediately exposes the Mq + Mi states in 
Ho and Hi as quasi-bound states with eigenvalue /in = 
(corresponding to instantaneous decay, 7 ^ 00). Does this 
exhaust all short-lived states? In view of degeneracy of the 
associated eigenvalue and the sub-unitarity of QFQ, this is 
by no means guaranteed. Indeed, Eqs. ^ naturally lead to a 
partial Schur decomposition [20], 



QFQ^UTU\ 



U={Po,Pu...,Pr^^,,P), 

Til Ti2\ (4) 

T22 



which reveals a much larger number of small eigenvalues. The 
unitary matrix U leaves the eigenvalues invariant. The struc- 
ture of the blocks of T is obtained by considering the opera- 
tion of QFQ on the columns in Pt. Equation ^ implies that 
Til is composed of sub-blocks that connect each subspace Ht 
to Ht-i (but not in the opposite direction). Hence, this ma- 
trix is upper triangular and all the diagonal elements vanish. 
It follows that the eigenvalue fin = indeed has an algebraic 
multiplicity 120] Mehr = Et<r^^, ^t- 



The difference between M^hr and Mq + Mi can be traced 
back to the break-down of the conditions for the conven- 
tional Weyl law known for closed systems |21], which as- 
sumes an orthogonal set of eigenvectors and translates this 
into a uniform covering of phase space. Perturbation theory 
shows that the small leakage out of the subspaces in Eqs. ^ 
lifts the degeneracy of eigenvalues /in ^ 0, but the result- 
ing MEhr eigenvectors are almost-linearly dependent, with 
their major component confined to the subspaces Hq and Hi. 
The conventional Weyl law would underestimate their num- 
ber as Mq + Mi ^ M{\Bo\ + \Bi\). The Schur decompo- 
sition carefully accounts for the admixture of the other sub- 
spaces Ht with t > 2, and provides an orthonormalized basis 
(Po, A, • • • , ^TEhr) fo^ the short-lived states. Equipped with 
such a basis, we now can formulate a new Weyl law in anal- 
ogy to the conventional case of a closed systems: The rank 
Mt of the characteristic projectors Vt can be estimated by the 
area covered by the regions Bt, Mt ~ M\Bt\. A direct con- 
sequence is the estimate MEhr ^ M\B\ for the number of 
short-lived states. 

Up to this point, our construction applies independent of 
the nature of the classical dynamics, which may be inte- 
grable, mixed, or chaotic. From now on we focus on clas- 
sically chaotic systems, which show a minimal amount of 
system- specific details provided the dynamics are 'sufficiently 
ergodic' (thwrIiA ;» 1). In parallel to the considerations 
in Ref. oO, a typical area (|A^,i|) ^ r^^^^^ex.p{-tX) 
shrinks exponentially with time t and equates to the Planck 
cell /leff = 1/M at TEhr = ^InA/". The number of in- 
dividual regions in Bt at given t is found by balancing the 
exponential shrinking of areas with the survival probability 
P{t) ~ T^gjj exp(— t/rdweii) for uniformly sampled starting 
points in phase space. The probability density P(|A|) to ini- 
tially reside in a region of area \A\ then takes the power law 
P{\A\) ^ A-^(|A|rdweii)"^"^/^^^^^^^ Collecting these re- 
sults, the region of quantum-to-classical correspondence cov- 
ers an area \B\ ^ 1 — exp(— TEhr/ Td well). We arrive at the 
estimate 



MEhr ~ M[l - exp(-rEhr/Tdwell) 



(5) 



of the number of instantaneous decay modes in a classically 
chaotic system. 

The complementary region B covers an area 



\B\ ^ exp(-rEhr/rdweii) = N' 



-l/Ardwell 



(6) 



In this region we can expect fully developed wave chaos, de- 
scribed by RMT with an effective matrix dimension 



M = M - MEhr ~ MA/--^/^^^--^ 



(7) 



Finally, with the RMT prediction of Ref. |22], the number of 
decay modes with 7^, smaller than a fixed value 7 < 1/rEhr 
is estimated as 



For fixed Tdweii we obtain n^ oc h^^, hence, a power-law 
dependence on /leff = 1/M with non-integer exponent d = 
(1/Ardweii) — 1- This is precisely of the form of the fractal 
Weyl laws recently put forward by Liu et al. |6], who for- 
mulated a trace formula for the long-lived quasi-bound states 
associated to the classical repeller of a chaotic open system 
(while we arrived at the result by constructing the short-lived 
states). 

We examine our predictions for a representative system 
with ballistic openings, the open kicked rotator also used in 
previous investigations |11, 12, 16] of the Ehrenfest time. In 
the classical limit, the phase space of this system is the torus 
(x^p) G [0, 1) [0, 1) with periodic boundary conditions. The 
phase-space density pt-\-i{x^p) = QJ^{pt{x^p)) evolves in 
discrete time-steps according to the standard map JT, speci- 
fied by 

K 

Xt-\-i = Xt -\- Pt -\ sin 27rxt (mod 1), 

47r 

K 

Pt+i =Pt^ -— (sin27rxt + sin27rxt+i) (mod 1), (9) 

47r 

which is followed by leakage in terms of a projection opera- 
tor Q that discards the density in a coordinate strip (x,p) G 
Ao = (xo, 3;o + l/rdweii)(0[0, 1). For kicking strength AT > 7, 
the standard map displays well-developed global chaos, char- 
acterized by a Lyapunov exponent A ^ ln(iir/2). We fix 
K = 7.5, xo = 0.2, Tdweii = 5. The quantum dynamics 
takes place in a Hilbert space of even integer dimension M, 
spanned by a basis of discretized position states atx = m/M, 
m = 0,l,...,M— 1. The time-evolution operator is 



r(- 



^r 



-+COS - 



(10) 



= ^[l-^d"Jell(l-e-")-']- 



(8) 



andQ = diag(]lArxAr, OatxAt, ll(M-2Ar)x(M-2Ar)). 

Husimi phase-space representations H{x^p) = 

Kx(^7P)IV^n)P of representative quasi-bound states for 
M = 160 are shown in Fig.[llb-f). The short-lived states in 
panels (b,c) have major component in the space Hi associated 
to the region Ai , and overlap in violation of the assumptions 
for the conventional Weyl law. Long-lived states (d,e) show 
random, delocalized wave patterns. Exceptions are a small 
number of long-lived trapped states, which are localized in 
the region of slowest escape. For M = 160 only one such 
state exists, which is shown in Fig.[ljf). 

Figure |2ja) shows the decay factors \pn\ = exp(— 7^/2), 
ordered by the decay rate 7^ (starting with the long-lived 
states). The data covers a range of inverse Planck's constants 
between M = 80 and M = 20480, and is presented as a 
function of the relative index n/M. In accordance with the in- 
creasing area covered by the region B of quantum-to-classical 
correspondence, the relative fraction of short-lived states in- 
creases with M. Figure |3b) shows the decay factors as a 
function of the renormalized relative index n/M, with M es- 
timated from Eq. Q and A = 1.32 approximated by the Lya- 
punov exponent of the closed system. The collapse onto one 
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FIG. 2: (Color online) Ordered decay factors |/in| = exp(— 7n/2) 
of the open kicked rotator with M = 2"^ • 80, m = 0, . . . , 8, and the 
other parameters as in Fig. [Ha), plotted as a function of the relative 
index n/M (panel a) or of the renormalized relative index n/M, 
where M is obtained from Eq. (panel b). The solid curve is the 
RMT result 1 22]. The dashed curve is this result with M fitted to the 
data. 



curve for all values of M confirms the predicted scaling Q of 
the number of short-lived states with the Ehrenfest time. 

Also shown in Fig.|2|is the result of RMT, obtained by sub- 
stituting the matrix F with a random matrix from the circular 
orthogonal ensemble. The almost perfect agreement found by 
fitting the effective matrix dimension M to the data [dashed 
curve in Fig. EJb)] demonstrates that effective RMT applies 
to a vast majority of the long-lived states. Since the number 
of states n^ is obtained by inverting this curves, the collapse 
of the data onto the effective RMT curve confirms the fractal 
Weyl law (Hi. 

In summary, we identified instantaneous decay modes (ex- 
ceptionally short-lived quasi-bound states) in open ballistic 
quantum systems, which capitalize on escape routes shorter 
than the Ehrenfest time. The large numbers in which these 
states emerge in the semiclassical limit is revealed only after a 
regularization of Weyl's law, which is required because of the 
almost-linear dependence of these states. For chaotic classical 
dynamics, the remaining long-lived states obey a RMT with 
an effective matrix dimension complying to a fractal Weyl 
law, from which one can conclude that such laws and the for- 
mation of instantaneous decay modes are intimately related. 

We gratefully acknowledge enlightening discussions with 
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